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ABSTRACT 

Using the Herschel Space Observatory's Heterodyne Instrument for the Far- 
Infrared (HIFI), we have observed para-chloronium (H 2 C1 + ) toward six sources in 
the Galaxy We detected interstellar chloronium absorption in foreground molec- 
ular clouds along the sight-lines to the bright submillimeter continuum sources 
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Sgr A (+50 km/s cloud) and W31C. Both the para-Hf C1+ and para-Hf Cl + iso- 
topologues were detected, through observations of their l n — o transitions at 
rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho- 
to-para ratio of 3, the observed optical depths imply that chloronium accounts 
for ~ 4 — 12% of chlorine nuclei in the gas phase. We detected interstellar chloro- 
nium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar 
photodissociation region and the Orion South condensation. For an assumed 
ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply 
total beam-averaged column densities of ~ 2 x 10 13 cm~ 2 and ~ 1.2 x 10 13 cm~ 2 , 
respectively, for chloronium in these two sources. We obtained upper limits on 
the para-H^Cl -1- line strengths toward H2 Peak 1 in the Orion Molecular cloud 
and toward the massive young star AFGL 2591. The chloronium abundances 
inferred in this study are typically at least a factor ~ 10 larger than the predic- 
tions of steady-state theoretical models for the chemistry of interstellar molecules 
containing chlorine. Several explanations for this discrepancy were investigated, 
but none has proven satisfactory, and thus the large observed abundances of 
chloronium remain puzzling. 

Subject headings: ISM: molecules - Submillimeter: ISM - Molecular processes 

1. Introduction 

With solar abundances estimated as 3.6 x 10~ 8 and 3.2 x 10~ 7 respectively relative 
to hydrogen (Asplund et al. 2009), fluorine and chlorine are only minor constituents of 
the Universe, with abundances 3-4 orders of magnitude smaller than those of carbon, 
nitrogen and oxygen. Nevertheless, these halogen elements show a strong tendency toward 
the formation of hydrides within the interstellar medium, with the surprising result that 
hydrogen fluoride (HF) and chloronium (H2C1 + ) can exhibit abundances that are comparable 
to those of CH and H 2 within diffuse molecular clouds (Neufeld et al. 2010a, Sonnentrucker 
et al. 2010, Lis et al. 2010). This tendency is a consequence of a unique characteristic of the 
thermochemistry of these elements: they are the only elements that can react exothermically 
with H-2 when in their predominant ionization stage in atomic clouds (i.e. atomic F, in the 
case of fluorine, and Cl + , in the case of chlorine). These reactions lead respectively to HF, 
which becomes the dominant reservoir of fluorine in diffuse molecular clouds (Neufeld et 
al. 2010; Sonnentrucker et al. 2010), and HC1 + , the tentative detection of which has been 
reported recently (De Luca et al. 2011). 

HF is relatively inert, being destroyed slowly via photodissociation and reaction with C + . 
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HC1 + , by contrast, can be destroyed by dissociative recombination, or by exothermic reaction 
with H 2 . The latter process leads to interstellar chloronium, discovered recently by Lis et al. 
(2010), who used the Herschel Space Observatory to perform absorption line spectroscopy of 
the In— Ooo transition of para-chloronium toward the bright submillimeter continuum sources 
NGC 66341 and Sgr B2 (S). Chloronium is an intermediary in the production of hydrogen 
chloride (Neufeld & Wolfire 2009, and references therein), the first chlorine-bearing molecule 
to be detected in the interstellar medium (Blake, Keene & Phillips 1985); the latter is formed 
as a result of dissociative recombination of H 2 C1 + , or - in dense clouds - following proton 
transfer from H2C1 + to CO. 

In this paper, we report four additional detections of chloronium obtained with Her- 
scheFs Heterodyne Instrument for the Far-Infrared (HIFI; de Graauw et al. 2010) in two 
separate Key Programs. We have observed the In — Ooo transition of para-H 2 Cl + , which lies 
less than 2 GHz from a molecular oxygen line that has been targeted in very deep integra- 
tions towards many sources in the HOP ("Herschel Oxygen Project"; Co-PI's, P. Goldsmith 
& R. Liseau) Open Time Key Program. Although studying H 2 C1 + was not the primary goal 
of HOP, we ensured that chloronium was observed simultaneously with 2 by a suitable 
choice of local oscillator setting, thus performing sensitive searches for chloronium towards 
all HOP sources. The HOP source list includes a wide diversity of molecular environments, 
and the new detections of H 2 C1 + reported here were obtained toward (1) the Orion Bar, 
a photodissociation region where intense radiation from the Trapezium stars in Orion is 
incident upon a dense molecular cloud; (2) Orion South, a dense condensation within the 
Orion Molecular ridge; and (3) the so-called "+50 km/s cloud" in the vicinity of Sgr A, a 
source of far-infrared continuum radiation - abutting the Sgr A East supernova remnant - 
with a sight-line that is intersected by multiple foreground molecular clouds located between 
the Sun and the Galactic Center. In addition, the In — Ooo transition of para-H 2 Cl + was 
targeted in observations of the massive star-forming region W31C, performed as part of the 
PRISMAS ("Probing Interstellar Molecules with Absorption Line Studies"; P.I., M. Gerin) 
Guaranteed Time Key Program. W31C (also known as G10.6-0.4) is a bright continuum 
source, located at an estimated distance of 4.8lo 8 kp c (Fish et al. 2003) within the so-called 
"-30 km/s" spiral arm, and provides a sight-line that intersects several foreground molecular 
clouds, the arrangement of which has been elucidated by Corbel & Eikenberry (2004; see 
their Figure 8). Because W31C shows a large continuum flux at infrared and radio wave- 
lengths, this sight-line has proven one of the most valuable in the Galaxy for the study of 
interstellar gas by absorption line spectroscopy. In addition to the new detections of chloro- 
nium, upper limits on the In — Ooo transition have been obtained from the HOP program 
toward H 2 Peak 1 in Orion, a source of shock-excited molecular hydrogen emission, and 
toward the young massive star AFGL 2591, proposed to be a strong far-ultraviolet emitter 
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that irradiates outflow cavity walls within a dense infalling envelope (Bruderer et al. 2010). 

The observations and data reduction methods are described in §2 below. In §3, we 
present the resultant spectra and report detections of H 2 C1 + toward four sources. In §4 and 
§5, the column densities and abundances inferred from these detections are discussed in the 
context of models for the chemistry of chlorine- containing interstellar molecules. 



2. Observations and data reduction 

All the observations reported here were obtained using the HIFI instrument in dual 
beam switch (DBS) mode, with the In — Ooo transition of para-chloronium in the lower 
sideband of the Band la receiver. Table 1 lists the observations that led to detection^. The 
telescope beam was centered at the coordinates given in Table 1, with the reference positions 
located at offsets of 3' on either side of each source. The beam size was 44" (half power beam 
width). The data were acquired using the Wide Band Spectrometer (WBS), which provides 
a spectral resolution of 1.1 MHz and a bandwidth of ~ 4 GHz. For each observation, the 
date and total duration are given in Table 1, together with the r.m.s. noise achieved in a 
1.1 MHz channel. As discussed in Neufeld et al. (2010a), the data were processed using 
the standard HIFI pipeline to Level 2, providing fully calibrated spectra with the intensities 
expressed as antenna temperature and the frequencies in the frame of the Local Standard 
of Rest (LSR). For the observations of W31C and Sgr A (+50 km/s cloud), we used three 
local oscillator (LO) frequencies, separated by a small offset, to confirm the assignment 
of any observed spectral feature to either the upper or lower sideband of the (double side 
band) HIFI receivers. The resultant spectra were coadded so as to recover the signal-to- 
noise ratio that would have been obtained at a single LO setting. Spectra obtained for the 
horizontal and vertical polarizations were found to be very similar in their appearance and 
noise characteristics and were likewise coadded. In the case of the Orion Bar and Orion S, 
the data were acquired in the frequency scan observing mode, which automatically performs 
the observations at multiple LO settings. Here, sideband deconvolution was performed using 
the doDeconvolution tool in the Herschel Interactive Processing Environment (HIPE). 

The relevant spectroscopy for chloronium is well established, thanks to direct laboratory 
measurements of its 180 - 500 GHz rotational spectrum (Araki et al. 2001). For the more 
abundant H^Cl" 1 " isotopologue, the strongest (F — 5/2 — 3/2) hyperfine component of the 
In — Oqo transition has a measured frequency of 485417.670 ± 0.015 MHz, with the F = 



1 Upper limits were obtained toward Orion H 2 Peak 1 (a = 5h35m 13.70 s, 6 = -5d 22' 09.00" J2000) and 
AFGL 2591 {a = 20h 29m 24.70 s, S = +40d 11' 19.00" J2000) 
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3/2 — 3/2 and F = 1/2 — 3/2 components lying at frequency offsets of -4.24 and +3.13 MHz, 
respectively, corresponding to Doppler velocities of +2.6 and -1.9 km/s. For H2 7 C1 + , the 
In - 00 F = 5/2 - 3/2 transition has a measured frequency of 484231.804 ± 0.026 MHz, 
and the F = 3/2 — 3/2 and F — 1/2 — 3/2 components lie at velocity offsets of of +2.1 and 
-1.7 km/s. The dipole moment of chloronium has been estimated as 1.89 Debye in ab initio 
calculations performed by Miiller (2008). 

3. Results 

Clear detections of the H2 5 C1 + In — Ooo transition were obtained in absorption toward 
Sgr A (+50 km/s cloud) and W31C, and in emission toward the Orion Bar and Orion South. 
In addition, the In — Ooo line of the Hij 7 Cl + isotopologue was detected toward Sgr A (+50 
km/s cloud) and W31C, but this transition was outside the spectral region covered by the 
Orion Bar and Orion S observations. 

Figure 1 shows the spectra of H2 5 C1 + (and, for two sources, HfCl + ) obtained toward 
the four sources Sgr A (+50 km/s cloud), W31C, the Orion Bar, and Orion South. Here, the 
intensity is plotted on the scale of antenna temperature, as a function of the LSR velocity 
for the strongest hyperfine component. For the Orion Bar and Orion S, Figure 1 shows 
the continuum-subtracted sideband-deconvolved spectra, whilst for Sgr A (+50 km/s cloud) 
and W31C, Figure 1 shows double sideband (DSB) spectra without continuum subtraction. 
In DSB spectra, the complete absorption of radiation in the observed spectral lines would 
reduce the antenna temperature to one-half its continuum value (given a sideband gain ratio 
of unity). The single sideband continuum antenna temperatures are listed in Table 1. For 
the Orion Bar and Orion South, the frequency-integrated line fluxes are also listed in units 
of K kms -1 , along with the line centroids and line widths. For the latter two parameters, 
two fitting methods were used: in the first method, the spectrum was fit as the sum of three 
Gaussians, of identical width, with the centroids separated in accord with the hyperfine 
splitting for the H2 5 C1 + transition, and with relative strengths appropriate to L.T.E.; in 
the second method, a single-component Gaussian fit was adopted, computed on the velocity 
scale for the strongest hyperfine component. The second method, of course, leads to line 
widths that are broader because of the effects of hyperfine splitting, and also results in a 
~ 0.8 kms -1 shift in the centroid velocity. In the Orion bar, the former method yields a 
line centroid (10.5 kms -1 ) that is in excellent accord with that observed for other spectral 
lines, although the line width (~ 3.4 kms -1 FWHM) is somewhat larger than those typically 
observed for optically-thin molecules along this sight-line (e.g. Parise et al. 2009), arguing, 
perhaps, in favor of a surface origin for H2C1 + . Toward Orion S, good agreement is obtained 
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with the line centroid (7.2 km s x ) and width (3.7 km s 1 ) inferred by Rodriguez-Franco et 
al. (2001) from observations of CN iV = 3 — 2. 

Toward W31C, several interloper lines are present in the H^Cl" 1 " spectrum. In particular, 
the S0 2 13 3)U - 12 2 ,i , S0 2 27i :27 - 26 , 26 and NS 2 n 1/2 J = 21/2 - 19/2 emission lines 
appear respectively near -25 +2, and +48kms _1 . The latter (NS 484.151 GHz) transition, 
in particular, affects our ability to determine the H 2 7 Cl + column density at LSR velocities 
above ~ 40kms~ 1 . However, this NS transition is a member of a lambda doublet, with 
another component appearing at 484.547 GHz (not shown in Figure 1). We have attempted 
to remove the interfering doublet member by subtracting a feature with the same strength 
and profile observed for the 484.547 GHz component; the result of this subtraction appears 
as the blue histogram in Figure 1, and has been used in our subsequent determination of the 
H^Cl" 1 " column density. 

Toward Orion H 2 Peak 1, the effective sensitivity was significantly compromised by the 
presence of several interfering spectral lines, the influence of which - at this position - is 
probably a consequence of the beam partially overlapping the nearby Orion hot core (a rich 
emission line source located at an offset of 27" from the beam center.) Based upon the typical 
residuals remaining after the removal of known transitions of dimethyl ether and 33 S0 2 , we 
obtained an upper limit of ~ 0.6 K km/s for the integrated H2 5 C1 + In — Ooo line intensity. 
Toward AFGL 2951, there is no evidence for emission or absorption at any Doppler velocity 
in a spectrum with a continuum antenna temperature of 57 mK and an r.m.s. noise of 4 mK 
(in a 1 km/s channel). 



4. Inferred column densities and abundances of chloronium 

4.1. Chloronium absorption observed toward W31C and Sgr A 

The Sgr A (+50 km/s cloud) and W31C spectra show H 2 C1 + absorption over a wide 
range of velocities, a behavior that is typical of other molecules observed in absorption toward 
these sources. The observed molecular absorption is believed to arise largely in foreground 
diffuse clouds, unassociated with the background sources of continuum radiation. 

In Figures 2 and 3, we present the observed flux for the two sources showing H^Cl 4 " 
and H| 7 C1 + in absorption, now normalized with respect to the continuum flux in a single 
sideband. Here, we assume a sideband gain ratio of unity (Roelfsema et al. 2012), and thus 
we plot the quantity 2T J 4/T J 4(cont) — 1. Using the procedure described by Neufeld et al. 
(2010b), we have deconvolved the hyperfine structure to obtain the spectra shown in red; 
these are the spectra that would have resulted if the hyperfine splitting were zero. For 
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comparison, the spectra (see figure captions for the sources of these spectra) of four other 
molecules detected with HIFI are shown: CH (the 536.761 GHz transition in green), HF (the 
1232.476 GHz transition in brown), para-H 2 (the 1113.343 GHz transition in magenta), and 
OH + (the 971.804 GHz transition in blue). In addition, the H I 21 cm spectrum is shown in 
dark green. Vertical offsets have been introduced for clarity. 

The spectra shown in Figures 2 and 3 indicate that all the absorbing species show broadly 
similar distributions in velocity space, but with HF, H 2 0, and CH - all believed to be good 
tracers of molecular hydrogen - showing the strongest mutual resemblance. Toward the Sgr 
A (+50 km/s cloud), chloronium is apparently distinctive among the molecules in exhibiting 
a lack of absorption at LSR velocities smaller than ~ — 70 km/s (and specifically a lack of 
absorption at ~ —130 km/s, which normally arises within the Expanding Molecular Ring 
(EMR) in the inner region of the Galaxy.) In this regard, the distribution of H 2 C1 + in velocity 
space is more similar to that of atomic hydrogen than that of any of the other molecules 
shown in Figure 3. Although, given the sensitivity of the observations, this behavior is only 
tentatively established, it argues in favor of theoretical predictions (§see 5.1 below) that the 
chloronium abundance is largest within clouds that are primarily atomic. 

In Table 2, the column densities of para-H 2 5 Cl + and para-H 2 7 Cl + are given for various 
velocity ranges, computed with the approximation that all para-chloronium molecules are in 
the ground Ooo state. As for other hydrides observed toward these sources, that approxima- 
tion is justified by the large spontaneous radiative rate for the observed transitional, the low 
density of the foreground clouds responsible for the absorption, and the weak continuum ra- 
diation field at the frequency of the observed transition and at the location of the absorbing 
material. Table 2 also presents the total H2C1 + column density (both isotopologues) for an 
assumed ortho-to-para ratio (OPR) of 3. For each velocity range, the H I column density is 
also given, determined from 21 cm absorption spectra (Fish et al. 2003 for W31C; Lang et 
al. 2010 for Sgr A +50 km/s cloud) for an assumed H I spin temperature of 100 K. Since 
the observed H2C1 + is expected to arise primarily in clouds of small H 2 fraction (see §5.1 
below), we also give the A r (H 2 Cl + )/A r (H) abundance ratios in Table 2. The mean value ob- 
tained for the to 50 km/s LSR velocity range in W31C, iV(H 2 Cl+)/iV(H I) = 1.15 x 10~ 8 , 
corresponds to ~ 12% of the gas-phase chlorine abundance within diffuse atomic clouds 
(1.0 x 10~ 7 ; Moomey, Federman & Sheffer 2011); that obtained for the -70 to -30 km/s LSR 
velocity range in Sgr A (+50 km/s cloud), A r (H 2 Cl + )/A^(HI) = 4.4 x 10~ 9 , corresponds to 



2 The Einstein A-coefficient for the observed para-H2 5 Cl + transition is 1.59 x 10~ 2 s -1 , for an assumed 
dipole moment of 1.89 D. This implies that the H2 and electron "critical densities" (at which the rate of 
collisional deexcitation is equal to this spontaneous radiative decay rate) are ~ 10 8 cm -3 and ~ 10 4 cm -3 
respectively, given the estimated collisional rate coefficients discussed in §4.2 
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~ 4% of the gas-phase chlorine abundance. For most velocity intervals, in both sources, the 
H2 5 C1 + /H2 7 C1 + abundance ratio is similar to the 35 C1 + / 37 C1 + isotopic ratio of 3 observed in 
the solar system (Lodders 2003). 

4.2. Chloronium emission observed toward the Orion Bar 

The absorption line observations discussed above provide robust estimates of the H 2 C1 + 
column densities, because almost all para-chloronium molecules are in the ground Ooo state 
(see §4.1 above), but the interpretation of the emission line observations is less straightfor- 
ward. At the densities of the H 2 C1 + emission region in the Orion Bar - n(H 2 ) ~ 3 x 10 4 cm -3 
and n(e) ~ 5 cm -3 (discussed further below) - the molecule will be subthermally-excited, 
with most para-H2Cl + molecules in the ground Ooo state, and thus the emission line in- 
tensity will depend crucially upon the rate coefficients for collisional excitation of H 2 C1 + . 
Because all excitations, to any excited rotational state, will be followed by a radiative cas- 
cade resulting in the emission of a In — Ooo hne photon, the relevant parameter is the total 
rate of collisional excitation out of the ground rotational state. As far as we are aware, 
no detailed calculations of the latter have been performed. For excitation by H2, we there- 
fore adopt the results obtained (Faure et al. 2007) for the isovalent para-H 2 molecule, 
which yield a total rate coefficient for excitation from Ooo that is well-approximated by 
gtot(H 2 ) = 2 x 10~ 10 T 2 a5 exp[- J E(l 11 )/A;T]cm 3 s _1 , where E(l n ) is the energy of the low- 
est excited state (In), T is the kinetic temperature, and T 2 = T/100 K. For the case 
of collisional excitation by electrons, we used the Coulomb-Born approximation to relate 
the spontaneous radiative rate for the In — Ooo transition to the rate of electron impact 
excitation. With the aid of equations 15, 12, and 11 in Neufeld & Dalgarno (1989), we 
thereby obtained an estimate for the rate coefficient for excitation from Ooo that can be fit 
by q tot (e) = 2.6 x 10~ 6 T 2 ^ a37 exp[- J E(ln)/A;T] cm 3 s" 1 , where E(l n ) = 23 K for the case of 
H 2 C1 + . We have also considered - and then rejected - the possibility that excited states of 
H 2 Ci + might be populated significantly when the molecule is formed. Even if every H 2 Ci + is 
formed in an excited state, this so-called "formation pumping" can be shown to be negligible 
by means of an argument that equates the rate of formation of H 2 C1 + to the rate of its 
destruction through dissociative recombination and proton transfer to neutral speciesjf] 



3 TheH 2 Cl+ destruction rate is k BR n e +k PT n(X), where fc DR , = 1.2 x 10" 7 (T/300K)-°- 85 cm 3 s" 1 (Geppert 
& Hamburg 2009, personal communication) is the rate coefficient for dissociative recombination of H2C1 + , 
kpT ~ 10 _9 cm 3 s _1 is a typical rate coefficient for proton transfer, and n(X) ~ few x 10 _4 n(H2) is the 
density of neutral species (e.g. CO) to which proton transfer might occur. The first term in this expression 
is smaller than the rate of electron-impact excitation, whilst the second term is much smaller than the rate 
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The molecular hydrogen density in the Orion Bar has been estimated previously by 
means of PDR (Photodissociation Region) modeling (e.g. Jansen et al. 1995). The consensus 
estimate, which we adopt here, implies a total density of hydrogen nuclei of toh ~ 6 x 
10 4 cm -3 , although previous studies have suggested the presence of denser clumps with a 
small filling factor (e.g. Lis & Schilke 2003); in the present work, we neglect the latter, 
because their contribution to the overall emission is expected to be small. Theoretical studies 
for the chemistry of chlorine-bearing molecules (e.g. Neufeld & Wolfire 2009, hereafter NW09) 
suggest that H 2 C1 + is most abundant near the cloud surface where carbon is fully ionized, 
hydrogen is largely atomic, and the gas temperature is ~ 500 K. In this region, the atomic 
hydrogen density is essentially the hydrogen nucleon density, n(HI) ~ nn, and we estimate 
the collisional excitation rate (per H 2 C1 + molecule) as -Rsurfaco = <Ztot(H2)n<H + 9tot(e)n e = 
3.8 x 10~ 5 s _1 . Here, we assumed that the rate of collisional excitation by H is equal to that 
for H 2 , and we adopted an electron abundance equal to the gas-phase abundance of carbon 
(~ 1.6 x 10~ 4 relative to H nuclei; Sofia et al. 2004). We also considered the conditions 
existing far from the cloud surface, in the cloud interior, where the fractional ionization 
is low, hydrogen is primarily molecular, and the gas temperature is ~ 30 K; if the H 2 C1 + 
originates primarily in such a region, the collisional excitation is -Ri nter ior = <?tot(H 2 )nH/2 = 
1.5 x 10~ 6 s~ 1 . Unfortunately, these estimates are quite uncertain, given the absence of 
reliable rate coefficients for H 2 C1 + . In particular, several shortcomings in the Coulomb- 
Born approximation have been discussed by Faure & Tennyson (2001), and there may be 
significant differences in the rates of excitation of H 2 C1 + and H 2 in collisions with H 2 . 

In the limit of subthermal excitation, the integrated brightness temperature of the 
In — Oqo line is related to the column density of para-H 2 Cl + by the expression 



where iV(p-H 2 Cl + ) is the beam averaged column density of para-H 2 Cl + . In deriving 
equation (1), we adopted the standard definition of Rayleigh- Jeans brightness tempera- 
ture, Tb = I u \ 2 /2k, where l v is the specific intensity; and we assumed that the line emis- 
sion is optically-thin and isotropic, with the frequency-integrated specific intensity given by 
/i^iV(p-H 2 Cl + )/47r. 

Given an integrated antenna temperature of 382 mK km/s for the In — Ooo transition of 
H 2 5 C1 + , and a main beam efficiency of 0.76 (Roelfsema et al. 2012), this expression yields an 
estimate of 7x 10 12 cm -2 for the beam-averaged column density of para-H 2 5 Cl + , if assumed to 




(1) 



of collisional excitation by H2 or H. 
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originate in a region close to the cloud surface (where R = -Rsurface)- Alternatively, if the para- 
H| 5 C1 + originates in the cloud interior (where R = -Rinterior); the estimated column density 
is even higher: 1.6 x 10 14 cm -2 . For an assumed ortho-to-para ratio of 3 for chloronium, 
corresponding to the equilibrium value expected at diffuse cloud temperatures, and for a 
H 2 5 C1 + /H 2 7 C1 + ratio of 3, these values correspond to total H2C1 + column densities of 2.9 x 
10 13 cm -2 and 7.0 x 10 14 cm -2 , respectively, for a cloud surface and cloud interior source of 
the observed H| 5 C1 + emission. These estimates both apply to the beam-averaged column 
density and represent two extreme cases that are expected to bracket the actual value. 



4.3. Chloronium emission observed toward the Orion South condensation 

In the Orion South molecular condensation, the gas density is believed to be higher than 
in the Orion Bar PDR (and UV radiation field to be smaller.) Mundy et al. (1988) estimated 
the virial mass of Orion S to be 80 M from interferometric CS J = 2 — 1 observations, and - 
assuming a thickness along the line-of-sight based on the observed angular size - obtained an 
estimate of 3 x 10 6 cm -3 for the H2 density, n(H 2 ) = wh/2. A better density estimate comes 
from fits to the J = 4 - 3, J = 10 - 9, J = 12 - 11 and J = 16 - 15 transitions of HC 3 N, 
observed by Bergin, Snell & Goldsmith (1996), who mapped the entire Orion ridge in these 
transitions with an angular resolution of about 50" . Assuming a gas temperature of 30 K, 
Bergin et al. obtained a density estimate n(H 2 ) > 1.2 x 10 6 cm -3 at the position of Orion 
S, with n(H 2 ) ~ 10 6 cm~ 3 at adjacent positions. For a higher assumed gas temperature 
~ 100 K, the density limit for Orion S would be lower by about a factor of 3. Finally, 
Rodriguez-Franco et al. (2001) combined observations of the N = 1 — 0, N = 2 — 1 and 
N = 3 — 2 lines of CN to estimate a density n(H 2 ) = 2 x 10 5 cm -3 for the region near Orion 
S, for an assumed gas temperature of 80 K. In our analysis of the H 2 C1 + emission detected 
from Orion S, we have adopted a density estimate ri(H 2 ) = 3 x 10 5 cm~ 3 , based on these 
previous values presented in the literature. 

The gas temperature in Orion South has been estimated from observations of two sym- 
metric top molecules: NH3 (Batrla et al. 1983) and CH3CN (Ziurys et al. 1990). Although 
Batrla et al. mapped the (1,1) and (2,2) inversion lines of NH3 with a resolution of 43" , only 
for a few select positions were observations of the higher inversion lines made and tempera- 
tures quoted. One such position was within 11" of Orion South (their position S6). These 
observations and those of Ziurys et al. are consistent with a gas temperature of about 100 
K for Orion S, a value that we adopt in deriving the H 2 C1 + column densitjQ. 



4 We note also that Mauersberger et al. (1986) detected the (J,K) = (7,7) metastable inversion line of 
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There being no clear evidence in Orion South for a PDR of the same prominence as that 
in the Orion Bar, we expect collisions with H 2 to dominate the excitation of chloronium. 
Given an integrated antenna temperature of 316 mK km/s, and the density and tempera- 
ture estimates discussed above, we derive a para-H^Cl" 1- column density of 3 x 10 12 cm -2 , 
corresponding to a total chloronium column density of 1.6 x 10 13 cm -2 for an assumed OPR 
of 3 and H 2 5 C1 + /H 2 7 C1 + ratio of 3. We estimate the beam-averaged H2 column density as 
4 x 10 23 cm~ 2 , based upon the 1.3 mm dust continuum map obtained by Mezger, Zylka & 
Wink (1990), which implies an N(R 2 C\ + )/N(R 2 ) abundance ratio of 4 x 10 -11 . This value 
corresponds to a fraction ~2x 10~ 4 of the gas-phase CI abundance determined for diffuse 
atomic clouds. By means of a similar analysis, we found that the upper limit on the H 2 C1 + 
line flux obtained toward Orion H 2 Peak 1 yields a limit of 5 x 10" 10 on the N(B. 2 C\ + ) /JV(H 2 ) 
abundance ratio at that position, a value that is a factor ~ 10 higher than the abundance 
inferred for the Orion South condensation. 



5. Comparison with model predictions 

The past four decades have seen the publication of many theoretical studies of the 
chemistry of chlorine-bearing molecules in the interstellar medium (Jura 1974; Dalgarno et 
al. 1974; van Dishoeck & Black 1986; Blake, Anicich & Huntress 1985; Schilke et al. 1995; 
Federman et al. 1995; Amin 1996; NW09). In diffuse molecular clouds, the chemistry is 
initiated by the photoionization of atomic chlorine by ultraviolet radiation in the narrow 
912 - 958 A wavelength range 

CI + hu -> C1+ + e, (Rl) 
which is followed by the exothermic reaction 

C1+ + H 2 HC1+ + H (R2) 

Depending upon the molecular hydrogen fraction, /(H 2 ) = 2ri(H 2 )/[2n(H 2 ) + n(H)], HC1 + 
is destroyed primarily via dissociative recombination 

HC1+ + e ->■ H + CI, (R3) 

or by reaction with H 2 to form H 2 C1 + 

HC1+ + H 2 -> H 2 C1+ + H. (RA) 



NH3, suggesting that in addition to the 100 K component, a much hotter component is also present with 
a temperature > 270 K, although based on the rotation diagram shown, less than 10% of the NH3 column 
density is in the high temperature component. 
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Chloronium ions are subsequently destroyed by dissociative recombination to form CI or 
HC1: 

H 2 C1 + + e -> CI + H 2 or 2H (ffia) 

H 2 C1 + + e -> HC1 + H (R5b) 

In dense molecular clouds, by contrast, where the interstellar UV radiation field is 
greatly attenuated, the reaction sequence is driven by cosmic ray ionization instead of pho- 
toionization. Here, cosmic ray ionization of H 2 leads to the reactive Hjj" ion (following proton 
transfer from to H 2 ), and reactions Rl and R2 are replaced by 

CI + H+ -> HC1+ + H 2 (RQ) 

as the primary route to HC1 + . 

For several key reactions in the above sequence, there are no reported measurements 
of the rate coefficients. In particular, the dissociative recombination (DR) rate for HC1 + 
an important parameter affecting the abundances of both HC1 + and H 2 C1 + - is unknown 
and urgently needed. In their recent study of chlorine chemistry in diffuse molecular clouds, 
NW09 adopted a rate coefficient for this process of 2 x 10~ 7 cm 3 s _1 , a value typical of DR rate 
coefficients for diatomic molecular ions; however, several hydride molecular ions, including 
OH + and HF + , have measured rate coefficients up to an order of magnitude smaller than 
this typical value. 



5.1. Diffuse cloud abundances toward W31C and Sgr A (+50 km/s cloud) 

To help interpret the H 2 C1 + observations of diffuse clouds along the sight-lines to Sgr 
A (+50 km/s cloud) and W31C, we have used the Meudon PDR code (le Petit et al. 2006) 
to obtain predictions for the H 2 C1 + abundance. The approach is rather similar to that 
adopted by NW09, but includes an improved treatment of the depth dependence of the 
photochemistry and of the self-shielding of molecular hydrogen from photodissociation. In 
particular, the photodissociation and photoionization rates for such species as H, H 2 , CI, 
HC1, HC1 + and CO have been computed by integrating the corresponding cross-sections 
over the radiation field intensity in an "exact" UV transfer approach where possible self- 
shielding, line overlap and continuum gas phase and dust absorption are taken into account. 
The contribution of the dust extinction is computed from the shape of the extinction curve 
whose analytic expression is given and discussed in Fitzpatrick & Massa (1999). The mean 
Galactic curve extinction has been used here as our standard assumption. The radiation 
field is assumed to impinge isotropically on both sides of a slab of fixed visual extinction 
Av,tot- The reaction network is the same as that adopted by NW09. 
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In Figure 4, we present the predicted fractions of gas-phase chlorine present in the 
HC1 + and H 2 C1 + ions, /(HC1 + ) and /(H 2 C1 + ), averaged along a sight-line through a diffuse 
molecular cloud of total visual extinction Ay jtot , together with the average H 2 fraction, 
/(H 2 ) = 2N(R 2 )/[2N(H 2 ) + iV(HI)]. Although the diffuse-cloud abundances of HC1+ and 
H 2 C1 + show very little dependence upon the cosmic ray ionization rate, we note that the 
results presented in Figure 4 were computed for an assumed primary cosmic ray ionization 
rate of 1.8 x 10~ 16 s _1 per H nucleus, a value corresponding to the "enhanced ionization rate" 
considered by NW09 and in agreement with recent estimates obtained from the analysis of 

(Indriolo et al. 2007) and OH + abundances (Neufeld et al. 2010b) in diffuse clouds. 
The different panels in Figure 4 show the results for different ratios of the gas density, 
nui to the radiation field, Xuv, the latter normalized with respect to the mean radiation 
field given by Draine (1978). Regardless of nn/xuv, the maximum value of /(H 2 C1 + ) is 
typically obtained in clouds with /(H 2 ) ~ 0.3. In Figure 5, we show (magenta diamonds) the 
ratio of the H 2 C1 + column density to that of atomic hydrogen, A r (H 2 Cl + )/A r (H I), together 
with the corresponding ratio for HC1 + (blue triangles). The maximum A r (H 2 Cl + )/A^(H I) 
ratio predicted for any parameters, 1 x 10~ 9 , is a factor of 4 and 12 below the average 
values measured toward Sgr A (+50 km/s cloud) and W31C, respectively. A very similar 
discrepancy has been reported in the case of HC1 + , the tentative detection of which has been 
presented recently by De Luca et al. (2011). Indeed, if we assume a Af(H 2 Cl + )/iV(HCi + ) 
ratio of 1, the value implied for W31C given the De Luca et al. detection, the maximum 
predicted A r (H 2 Cl + )/A^(H I) ratio is 5 x 10~ 10 , a factor of ~ 20 below the measured value. 

We have investigated several explanations for the discrepancies described above, none 
of which appears to be entirely satisfactory. First, we have examined the effect of reducing 
the rate coefficient assumed for DR of HC1+ to 2 x 10~ 8 cm 3 s~ 1 (T/300 K)^ 1 / 2 ; this value is 
10 times smaller than that adopted by NW09 as the typical DR of diatomic molecular ions, 
but equal to the experimental value obtained by Djuric et al. (2001) for HF + . The result of 
this change is quite modest, however: for the cloud model that yields the maximum H 2 C1 + 
abundance, with Ay tto t = 0.3 and nn/xuv = 32, the predicted HC1 + and H 2 C1 + abundances 
increase by factors of 2.7 and 1.3 respectively. Second, we investigated the effect of increasing 
the rate coefficient assumed for charge transfer between CI and H + . That reaction was 
assumed by Blake et al. (1986) to have a rate coefficient of 10 -9 cm 3 s -1 and was identified as 
an important route to the formation of Cl + , but subsequently found to have a considerably 
smaller rate coefficient (by factors of ~ 20 and ~ 250 at 300 K and 100 K, respectively) in 
calculations performed by Pradhan & Dalgarno (1994). For the cloud model with Ay t tot = 0.3 
and uh/xuv — 32, the predicted HC1 + and H 2 C1 + abundances only increase by 7% and 6%, 
respectively, if a rate coefficient of 10~ 9 cm 3 s _1 is assumed in place of that computed by 
Pradhan & Dalgarno (1994). Third, we have examined the effect of adopting an interstellar 
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radiation field with a spectral shape corresponding to a 30,000 K blackbody in place of that 
assumed by Draine (1978) (but with the absolute intensity constrained to be equal at 1000 
A .) At photon energies above the ionization potential of atomic chlorine (12.97 eV), the 
former drops less rapidly than the latter, leading to an enhanced photoionization rate for 
CI. However, even for this harder (and arguably extreme) radiation field, the resultant HC1 + 
and H 2 C1 + abundances are only increased by factors of 1.4 and 1.3 respectively. Fourth, we 
have investigated the dependence of the model predictions upon the wavelength dependence 
of the dust opacity, assumed in our standard model to follow the average Galactic extinction 
law of Fitzpatrick& Massa (1999). We have considered a case in which the dust opacity is 
assumed to rise less rapidly as the wavelength decreases, following the extinction curve of 
HD 294264 (Fitzpatrick & Massa 2005). Here, the peak H2C1 + abundances are achieved at 
somewhat larger Av,tot, but the effect on the peak abundances are small: the peak predicted 
HC1 + and H 2 C1 + abundances are only increased by 3% and 12% respectively. 

We note two additional systematic uncertainties that might help account for the dis- 
crepancy between theory and observation. First, the chloronium column densities derived 
from the absorption line observations are inversely proportional to the square of the as- 
sumed dipole moment, for which we adopted a value of 1.89 D, the result of unpublished 
calculations by Miiller. This value is already fairly large, and thus an unusually large dipole 
moment (larger than ~ 4 — 5 D) would be needed to bring the theory and observation into 
acceptable agreement. Second, the chloronium column densities predicted by theory are 
roughly inversely proportional to the assumed DR rate for H 2 C1 + . The value we adopted 
was that obtained in recent laboratory measurements ( Geppert & Hamberg 2009, personal 
communication) performed with the use of an ion storage ring, but there is always the 
possibility that some degree of H2C1 + rotational excitation in the laboratory environment 
may lead to results that differ from those applicable in the interstellar environment (where 
most molecules are in the ground rotational state); in this regard, we note that Petrignani 
et al. 2011 have recently argued that studies of dissociative recombination in ion storage 
rings may typically involve ions with higher rotational temperatures than had previously 
been appreciated. 

Finally, we note that the theoretical models employed to date obtain steady-state abun- 
dances for the chlorine-bearing molecules under investigation. Because the timescale for H 2 
formation is relatively long - of order 10 (100 cm~ 3 /nH) Myr - time-dependent effects may 
be important as molecular clouds either form or disperse. In the case of molecular cloud 
formation, the H 2 abundance lags relative to its equilibrium value, which is unfavorable for 
the formation of chloronium; but the case of cloud dispersal may be more favorable for the 
production of enhanced H 2 C1 + abundances and is worthy of future investigation. 
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5.2. The column density in the Orion Bar PDR 

The chloronium column density predicted by NW09 (their Figure 5) for conditions ap- 
propriate to the Orion Bar PDR - nu ~ 6 x 10 4 cm -2 and xuv ~ 3 x 10 4 (Herrmann et al. 
1997) - is iV(H 2 Cl + ) ~ 3 x 10 11 cm -2 , a factor ~ 100 smaller than the observed value (see 
§4.2, for the case of a surface source of H 2 C1 + ). These two column densities are not directly 
comparable, however, because the theoretical prediction is of the perpendicular column den- 
sity of chloronium, whereas the Orion Bar PDR is viewed nearly edge-on. Nevertheless, the 
geometric enhancement in the measured column density is expected to be a factor ~ 4 (e.g. 
Neufeld et al. 2006), and thus the theory still appears to underpredict the observed column 
density by more than an order of magnitude. Clearly, the uncertainties here are considerably 
larger than those involved in the comparison of theory and observation in diffuse clouds, but 
nevertheless it appears that the discrepancy is similar in size (and direction). 

5.3. The abundance in the Orion South condensation 

Whereas the chemistry leading to chloronium in diffuse clouds is driven by photoion- 
ization of CI, the chemical network in dense clouds is initiated by cosmic ray ionization. In 
Figure 6, we show theoretical predictions for the H^Cl" 1 " abundance in the shielded interior of 
a dense cloud, shown here relative to that of gas-phase chlorine. These predictions are based 
upon an assumed cosmic-ray ionization rate of 1.8 x 10 _17 s _1 (primary ionization rate per 
H nucleus); this is the "standard value" adopted by NW09, believed to be most appropriate 
for dense molecular clouds. Once again, the observed H 2 C1 + abundance of ~ 1.5 x 10 -4 
relative to the diffuse cloud gas-phase chlorine abundance (§4.3) is more than an order of 
magnitude greater than the model prediction. The discrepancy becomes even greater if the 
chlorine depletion in the Orion South condensation is larger than that in the diffuse ISM, as 
suggested by observations of HC1 by Schilke et al. (1997). 

6. Summary 

1. We have detected interstellar chloronium absorption in foreground molecular clouds 
along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) 
and W31C. Both para-H^Cl 4 " and para-H| 7 Cl + were detected, in absorption, in the In — o 
ground-state transition. For an assumed ortho-to-para ratio of 3, the observed optical depths 
imply that chloronium accounts for ~ 4 — 12% of chlorine nuclei in the gas phase. 

2. We have detected interstellar chloronium emission from two positions in the Orion 
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Molecular Cloud I: the Orion Bar photodissociation region and the Orion South condensa- 
tion. For an assumed ortho-to-para ratio of 3, the observed emission line fluxes imply total 
beam-averaged column densities of ~ 2.9 x 10 13 cm -2 and ~ 1.6 x 10 13 cm -2 , respectively, 
for chloronium in these two sources. 

3. The chloronium abundances inferred in this study are typically at least a factor ~ 10 
larger than the predictions of theoretical models for the chemistry of interstellar molecules 
containing chlorine. Several explanations for this discrepancy were investigated, but none has 
proven satisfactory, and thus the large observed abundances of chloronium remain puzzling. 
Additional laboratory and/or theoretical investigations of the dipole moment of H^Cl" 1 ", its 
dissociative recombination rate (for rotationally-cold ions), and its rates of collisional exci- 
tation by H 2 and e are urgently needed, as are time-dependent astrochemical models. 

HIFI has been designed and built by a consortium of institutes and university depart- 
ments from across Europe, Canada and the United States under the leadership of SRON 
Netherlands Institute for Space Research, Groningen, The Netherlands and with major con- 
tributions from Germany, France and the US. Consortium members are: Canada: CSA, 
U. Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, MPS; 
Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Osservatorio Astrofisico di Arcetri- INAF; 
Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astronomico Na- 
cional (IGN), Centro de Astrobiologa (CSIC-INTA). Sweden: Chalmers University of Tech- 
nology - MC2, RSS & GARD; Onsala Space Observatory; Swedish National Space Board, 
Stockholm University - Stockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: 
Caltech, JPL, NHSC. 

Support for this work was provided by NASA through an award issued by JPL/Caltech. 
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Table 1. Summary of observations 



Source 




Orion Bar 


Orion S 


W31C 


Sgr A 
(+50 km/s cloud) 


R.A. (J2000) 




5h 35m 20.60s 


5h 35m 13.50s 


18h 10m 28.70s 


17h 45m 51.70s 


Declination (J2000) 




-5d 25' 14.00" 


-5d 24' 8.00" 


-19d 55' 50.00" 


-28d 59' 09.00" 


Date 




2010-03-01 


2010-03-01 


2011-04-08 


2010-09-16 


Total duration (s) 




4176 


4176 


588° 


7848" 


Continuum T A (SSB in mK) 


66 


800 


382 


115 


Integrated line strength'' 


(mK km/s) 


328 


316 






Line centroid b (km/s) 




10.5/11.3 C 


7.8/8.5 c 






Line FWHM' (km/s) 




3.4/5. l c 


4.0/5. 7 C 






R.m.s noise (mK) 




4.5 


4.5 


9.7 


1.4 


Mode 




Frequency scan 


Frequency scan 


Pointed 


Pointed 


Key program 




HOP 


HOP 


PRISMAS 


HOP 



a Divided into three observations of equal duration with different LO settings (see text) 
b Emission lines only 

c The first number is for a three-component Gaussian fit, with the velocity components separated in accord with the 
hypcrfinc splitting and the relative strengths appropriate to L.T.E.; the second number is for a single- component Gaussian 
fit, computed on the velocity scale for the strongest hyperfine component. 



Table 2. Chloronium column densities and abundances from absorption line observations 



Source 


ulsr range 
(km/s) 


A r (p-H| 5 Cl+) 
10 12 cm- 2 


A r (p-H| 7 Cl+) 
10 12 cm- 2 


JV(H 2 C1+) 
10 12 cm- 2 


35 C1/ 37 C1 


7V H 
10 21 cm" 2 


AT(H 2 C1+)/7Vh 
lO" 9 


W31C 


[10, 20] 


3.9 


1.6 


22 


2.5 


2.2 


10.0 




[20, 30] 


5.9 


3.0 


36 


2.0 


2.9 


12.1 




[30, 40] 


8.5 


3.2 


17 


2.6 


3.5 


13.4 




[40, 50] 


2.8 


0.9 


15 


3.3 


1.7 


8.0 


Total: 


[10, 50] 


21.1 


8+ 


118 


2.5 


10 


11.5 


Sgr A 


[-70, -30] 


2.3 


0.9 


13 


2.5 


2.9 


4.4 




[-30, +20] 


5.5 


1.3 


27 


4.2 


> 14 


< 2.0 


Total: 


[-70, +20] 


7.8 


2.2 


40 


3.5 


> 17 


< 2.4 



a The p-H| 7 Cl + absorption spectrum in this velocity range is significantly affected by an interloper absorption feature 
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Fig. 1. — Spectra of the para-ln — Ooo transition (and, for two sources, the analogous 
transition of para-H| 7 Cl + ) obtained toward the four sources Sgr A (+50 km/s cloud), W31C, 
the Orion Bar, and Orion South. Here, the intensity is plotted on the scale of antenna 
temperature, as a function of the LSR velocity for the strongest hyperfme components (at 
rest frequencies of 485.417 GHz for Hf C1+ and 484.232 GHz for HfCl+ ) The spectra 
shown for Sgr A (+50 km/s cloud) and W31C are double sideband spectra, while those for 
the Orion Bar and Orion South are sideband deconvolved. 
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W31 C 

t | i i i | i i i | i i r 




LSR velocity (km/s) 

Fig. 2. — H| 5 C1 + and H2 7 C1 + absorption spectra obtained toward W31C, normalized with 
respect to the continuum flux in a single sideband. The black curves show the observed 
spectra, and the red curves show the result of deconvolving the hyperfine structure; these 
are the spectra that would have resulted if the hyperfine splitting were zero. For compar- 
ison, the spectra of four other molecules detected with HIFI are shown: CH (the 536.761 
GHz transition in green from Gerin et al. 2010b, after deconvolution of the hyperfine struc- 
ture), HF (the 1232.476 GHz transition in brown, from Neufeld et al. 2010a), para-H20 
(the 1113.343 GHz transition in magenta, from Neufeld et al. 2010a), and OH+ (the 971.804 
GHz transition, from Gerin et al. 2010a, both with (blue) and without (black) hyperfine 
deconvolution). In addition, the H I 21 cm spectrum is shown in dark green (from Fish et 
al. 2003). Vertical offsets have been introduced for clarity. 
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Sgr A ( + 50 km/s cloud) 




LSR velocity (km/s) 

Fig. 3. — Same as Figure 2, except for Sgr. A (+50 km/s cloud). Here, the H I 21 cm 
spectrum is from Lang et al. (2010), the OH + spectrum is from Neufeld et al. (2011), and 
the HF and H2O spectra are from Sonnentrucker et al. (2011). The HF and H2O spectra 
were obtained at a position 44" West and 20" South of that observed here. 
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Fig. 4. — Predicted fractions of gas-phase chlorine present in the HC1 + (blue trian- 
gles) and H 2 C1 + (magenta diamonds) ions, averaged along a sight-line through a diffuse 
molecular cloud of total visual extinction Ay .tot, together with the average H2 fraction, 
/(H 2 ) = 2iV(H 2 )/[2iV(H 2 ) + iV(HI)] (black). The different panels show the results for dif- 
ferent ratios of the gas density, %, to the radiation field, XuVi the latter normalized with 
respect to the mean radiation field given by Draine (1978). 
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lo gio( A v,tot/ ma g) | ogi (A v ,tot/ ma g) 

Fig. 5.— Predicted iV(HCl + )/iV(H I) (blue triangles) and iV(H 2 Cl + )/iV(H I) (magenta dia- 
monds) ratios, averaged along a sight-line through a diffuse molecular cloud of total visual 
extinction Ay, tat- The different panels show the results for different ratios of the gas density, 
riH) to the radiation field, Xuv, the latter normalized with respect to the mean radiation 
field given by Draine (1978). 
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3 4 5 6 7 

log 10 (H 2 density in cm -3 ) 

Fig. 6. — Theoretical predictions (from NW09) for the dense cloud H 2 C1 + abundance, shown 
here relative to that of gas-phase chlorine, for an assumed cosmic-ray ionization rate of 
1.8 x 10 s -1 (primary ionization rate per H nucleus) 



